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GRAPHICAL  USER  INTERFACES 
FOR  DIGITAL-IMAGE  FREQUENCY  ANALYSIS 


Introduction 

Frequency  analysis  plays  an  important  role  in  the  area  of  visual  research.  The  application 
described  in  this  technical  memorandum  provides  a  graphical  user  interface  that  simplifies 
loading,  analyzing,  and  storing  the  results  of  common  frequency  analysis  tasks.  The 
application  allows  the  user  to  load  image  files  and  compute: 

•  A  2-D  discrete  Fourier  transform  (DFT)  of  the  image 

•  A  polar  transform  (orientation  and  radial  spatial  frequency)  of  the  DFT 

•  A  1-D  modulation  versus  spatial  frequency  plot  produced  by  averaging  over 
specified  orientations  of  the  polar  transform 

The  analysis  outlined  above  can  also  be  perfonned  on  a  subsection  of  the  input  image.  The 
subsection  of  the  input  image  can  also  be  multiplied  by  a  Gaussian  window  to  reduce  edge 
effects. 


Installing  and  running  the  application 

This  program  is  designed  to  run  on  a  PC  (Windows)  with  or  without  Matlab. 

Running  the  application  with  Matlab 

To  run  the  application  with  Matlab,  simply  copy  the  entire  folder,  m-codes,  (see  attached 
CD)  to  any  desired  location  on  an  end-user  machine.  The  user  can  either  change  the  current 
directory  to  m-codes  or  use  the  setpath  command  in  Matlab  to  reference  the  application 
directory.  Type  “guifft”  on  the  command  line  in  Matlab  to  start  the  application. 

Running  the  application  without  Matlab 

To  run  the  application  without  Matlab,  it  is  necessary  to  first  install  Matlab  runtime 
software.  The  runtime  modules  are  installed  by  running  MCRInstaller.exe  which  can  be 
found  on  the  attached  CD.  After  installing  the  runtime  modules,  the  application  is  started 
by  double-clicking  on  guifft.exe  located  inside  the  folder  labeled  StandAloneExecutable. 
Running  the  application  using  Matlab  is  primarily  of  interest  to  Matlab  programmers  who 
wish  to  examine  or  modify  the  m-code.  In  addition,  running  the  application  using  Matlab 
does  not  require  installing  the  runtime  modules. 
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Analyzing  a  Full  Image 

After  starting  the  application,  a  window  will  be  displayed  on  the  screen  as  shown  Figure  1. 
There  is  a  User  Control  Panel  (UCP)  near  the  bottom  of  main  window  that  allows  the  user 
to  select  the  input  image  and  initiate  the  frequency  analysis.  The  left  section  of  the  UCP 
provides  a  preview  of  the  input  image.  The  size  of  the  image  in  pixels  is  displayed  below 
the  figure.  The  user  can  instruct  the  program  to  simulate  monitor  gamma  by  using  the 
check-box  and  the  associated  text  box  positioned  below  the  preview  image.  If  the  Gamma 
box  is  checked,  the  numerical  values  of  each  image  pixel  are  modified  according  to  the 
following  formula: 

j  _  j  Gamma 

1  out  1  In  , 

where  Gamma  is  the  value  entered  into  the  gamma  text  box.  If  the  gamma  checkbox  is 
unchecked,  the  values  in  the  input  files  will  be  used  for  the  frequency  analysis. 

The  center  section  of  the  UCP  lists  the  content  of  the  current  directory.  Clicking  on 
will  move  the  current  directory  up  one  level  and  clicking  a  folder  icon  will  make  that 
folder  the  current  directory.  The  user  can  filter  the  type  of  image  files  displayed  by  using 
the  pull-down  menu  below  the  directory  listing.  The  user  can  preview  an  image  by 
selecting  (single-click)  the  filename  in  the  folder  list.  The  filename,  “flowers.bmp,”  is 
selected  for  preview  in  Figure  1 .  The  program  will  only  process  files  with  the  formats 
shown  in  Table  1.  Selecting  other  files  will  cause  “Invalid  Type”  to  be  displayed  on  the 
preview  window. 

The  Perform  FFT  button  initiates  the  data  analysis.  Upon  completion,  a  message  box 
will  be  displayed  indicating  that  the  analysis  is  complete  along  with  the  value  of  the 
gamma  used  (Figure  2).  Select  the  OK  button  to  dismiss  the  message  box. 


Format 

Extension 

Comment 

Bitmapped  Image 

.bmp 

Windows  bitmapped  image.  It  must  be  a  24-bit  depth 
image. 

Flexible  Image 
Transport  System 

.fit 

Fonnat  used  for  the  transport,  analysis,  and  archival 
storage  of  scientific  data  sets  in  the  astronomical 
community.  [Used  by  CCD  cameras  supplied  by  the 
Santa  Barbara  Instrument  Group.] 

Luminance  Image 

.lum 

A  simple  binary  format  with  the  first  two  16  bit  words 
being  the  X  and  Y  dimensions  of  the  luminance  array, 
followed  by  the  array  of  single-precision  (32  bit)  floating 
point  values.  [Used  by  the  Lumetrix  IQCam  CCD 
camera.] 

Text  image 

.txt 

A  tab  separated  text  file.  A  row  of  the  text  file  should 
equal  a  row  of  the  image. 

Table  1:  Valid  Image  Formats 
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Figure  1:  Perform  FFT  on  a  Full  Image 
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Job  Done!!!!!! 


on® 


The  DFT  is  performed  using  the  FFT2  function  supplied  by  Matlab: 


F(p, q ) 


M  - 1 


I 

m  =  0 


N-l 

v  -j(2K/M)pm  -j(2n/N)qn 
>  f(m,n)e  e 


n  =  0 


P  =  0,1 . M- 1 

q  =  0,  1,  1 


where  p  and  q  are  the  indices  of  the  DFT  coefficients  and  m  and  n  are  the  indices  of  the 
input  image.  The  program  displays  the  (Logio)  magnitude  of  the  complex  valued  Fourier 
coefficients, 


M(p,q)  =  ^Freal{p,q)2  +Fimag(p,q)2 


Results  of  the  analysis  are  shown  in  four  figures  that  appear  above  the  UCP  (see  Figure  1). 


1)  The  input  image  in  the  upper  left  of  the  figure, 

2)  The  magnitude  of  the  2-D  discrete  Fourier  transform  (DFT)  of  the  image  in  the 

lower  left  of  the  figure, 

3)  A  polar  transformation  of  the  magnitude  of  the  DFT: 


theta  =  arctan 


( V 


frequency  =  yj  p1  +  q2 

where  theta  is  the  angle  between  a  radial  line  segment  starting  from  the  DC  point 
(center  point)  of  DFT  and  a  line  segment  starting  at  DC  and  continuing  along  the 
positive  horizontal  axis.  If  a  point  on  a  line  does  not  fall  on  a  frequency  bin,  a 
bilinear  interpolation  is  used  to  calculate  the  value  for  the  point.  Orientations 
were  sampled  in  one  degree  increments  and  spatial  frequency  was  sampled  at  the 
bin  width  of  the  DFT.  Note  that  the  first  value  of  radial  frequency  will  be  DC 
value. 

4)  A  1-D  modulation  verses  spatial  frequency  plot  produced  by  averaging,  at  a  given 
spatial  frequency,  over  all  orientations  of  the  polar  transfonn  data.  The  data  in  the 
plot  are  normalized  to  the  DC  value  of  the  DFT. 
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Before  clicking  the  Perform  FFT  button,  the  user  should  decide  whether  the  results  are  to 
be  saved.  The  four  figures  discussed  above  will  be  saved  using  the  Portable  Network 
Graphics  (PNG)  format  if  the  relevant  checkbox,  in  the  right  portion  of  the  UCP,  is 
checked.  Likewise,  the  numerical  values  produced  by  orientation  averaging  will  be  saved 
as  a  Microsoft  Excel  file  if  the  relevant  checkbox  is  checked. 

The  final  processing  option  available  in  the  UCP  is  labeled  Define  Subsection  of  a  Image. 
We  will  next  describe  the  details  of  this  process. 

Analyzing  a  Selected  Portion  of  an  Image 

The  user  can  select  a  portion  of  the  input  image  to  be  processed  by  clicking  on  the  Define 
Subsection  of  a  Image  button  in  the  UCP  (again,  see  Figure  1).  When  this  is  done,  a 
window  entitled  “ Define  Rectangular  Subsection  for  FFT  Analysis  ”  is  displayed  (see 
Figure  3). 


Figure  3:  Define  Rectangular  Subsection  to  be  cut  out 
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The  user  can  use  the  mouse  to  move  the  curser  to  a  desired  location  and  click  to  specify 
that  location  as  a  center  point  for  the  rectangular  portion  of  the  image  to  be  analyzed.  The 
center  point  can  also  be  specified  by  manually  entering  coordinates  at  the  lower  right  of  the 
window.  The  size  of  the  portion  of  the  image  to  be  analyzed  must  be  specified  manually 
using  the  text  boxes  at  the  lower  left  of  the  windows.  A  yellow  circle  displayed  on  the 
image  indicates  the  previous  center  point  selected,  and  a  red  circle  indicates  the  current 
center  point,  ato  change  the  location  of  the  current  center  point,  the  user  should  press  the 
button  labeled  Re-select  Center  Point  and  use  either  the  cursor  or  manual  entry  to  select  a 
new  location.  If  the  center  point  is  chosen  near  the  edge  or  corner  of  the  image  and  the 
size  of  the  rectangular  area  is  set  too  large,  a  window  containing  an  error  message  will 
appear  (see  Figure  4).  In  that  case,  it  will  be  necessary  to  either  choose  another  center 
point  or  choose  a  different  size  for  the  rectangular  area  so  that  rectangular  area  will  fit 
inside  of  image. 


Figure  4:  Image  out-of-range  warning  box 

After  the  desired  center  point  and  size  are  chosen,  the  rectangular  area  can  be  cut 
out  by  pressing  the  button  labeled  CUT.  A  GUI  labeled  Perform  FFT  on  Subsection  of  a 
Image  will  appear  and  will  include  the  image  portion  selected  (see  Figure  5).  There  are 
three  checkboxes  on  this  GUI.  If  the  box  labeled  Gaussian  Window  is  checked,  a  Gaussian 
of  the  form: 


1  (x-x0  f 

<N 

O 

1 

l2  4  J 

i2  4  J 

will  be  multiplied  with  the  image  prior  to  calculation  of  the  DFT.  The  standard  deviations 
along  the  x  and  y  dimensions  (sx  ,  sy)  are  chosen  so  that  the  pixel  values  at  the  edge  of  the 
image  are  approximately  zero.  This  procedure,  called  data  windowing,  will  have  the  effect 
of  reducing  spectral  leakage  which  is  the  technical  tenn  used  to  describe  the  loss  of 
amplitude  of  a  given  frequency  to  other  frequency  bins  in  the  DFT. 

As  was  the  case  with  the  analysis  perfonned  on  the  full  image,  the  appropriate  boxes  must 
be  checked  in  order  to  save  the  image  (in  this  case  both  with  and  without  the  Gaussian 
window),  the  FFT  output,  the  polar  transform  of  the  FFT  output,  and  the  spectrum  obtained 
by  orientation  averaging.  The  results  will  be  saved  with  names  that  are  generated  by  the 
application,  and  which  indicate  the  center  point  location,  the  size  of  the  cut  image,  and 
whether  the  process  is  done  with  Gaussian  window. 
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Figure  5:  Perform  FFT  on  Subsection  of  a  Image 
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Key  Method  Used 

1 .  Gamma: 


newlum  =  oldlum . Agamma 


2 .  Orientation  Averaging: 

nPAD=min (size (lum) ) ;  %lum  is  the  input  image  data 
z=evennumber (nPAD) ; 
if  z.  1 

center  =  nPAD/2+1; 
step=nPAD/ 4 ; 

Rho=nPAD/2;  %length  of  radius 

else 

center= (nPAD+1 ) / 2 ; 
step=round( (nPAD+1) /4) ; 

Rho=center;  %length  of  radius 

end 

%--  Polar  transformation  and  interpolation  (i  is  in  degree) 
for  i  =  1:180  %i  is  in  degree 
Theta  =  i*pi/180; 

[x,  y] =pol2cart (Theta, Rho) ; 

%%  [x, y, value] =improfile (I, [vXstart, vXend] , [vYstart, vYend] , N 
equally  spaced  points Interpolation  Method'); 
if  i  <  90 

[cx, cy, c] =improfile (fft2_Im, [center, center+x-1] , [center, center- 
y] ,Rho, 'bilinear' ) ; 
elseif  i==90 

[cx, cy, c] =improfile (fft2_Im, [center, center] , [center, center- 
y] ,Rho, 'bilinear' ) ; 
else 

[cx, cy, c] =improf ile ( f f t2_Im, [center, center+x] , [center, center- 
y] ,Rho, 'bilinear' ) ; 
end 

%  Matrix  of  Theta  vs.  Radius 
radius_line (i, : ) =c;  %  bin  1  is  the  DC 

end 

3.  Gaussian: 

s=size (handles . lumcrop) ; 

hi =f special ( ' gaussian ',[s(l)  l],s(l)/3); 

h2=f special ( ' gaussian ',[1  s(2)],s(2)/3); 
h=hl*h2 ; 

h=h . /max (max (h) ) ;  %Normalized 
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Contents  of  the  Attached  CD 

The  CD’s  contents  are: 

1 .  an  electronic  copy  of  this  document  in  PDF  format, 

2.  a  folder  that  contains  Matlab  m-codes  of  this  program, 

3.  a  folder  that  contains  all  the  stand-alone  executable  files, 

4.  a  folder  that  contains  the  files  flower.bmp,  radial_sinewave.bmp,  and  the  testimage 
in  .lum,  .txt  and  .fit  formats, 

5.  MCRInstaller.exe. 
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